home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
CD ROM Paradise Collection 4
/
CD ROM Paradise Collection 4 1995 Nov.iso
/
program
/
nrpas13.zip
/
MIDPNT.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
689b
|
36 lines
PROGRAM d4r5(input,output);
(* driver for routine MIDPNT *)
CONST
nmax=10;
VAR
a,b,s : real;
i,glit : integer;
FUNCTION func(x: real): real;
(* Function for testing integration *)
BEGIN
func := 1.0/sqrt(x)
END;
FUNCTION fint(x: real): real;
(* Integral of 'func' *)
BEGIN
fint := 2.0*sqrt(x)
END;
(*$I MODFILE.PAS *)
(*$I MIDPNT.PAS *)
BEGIN
a := 0.0;
b := 1.0;
writeln;
writeln ('Integral of func computed with MIDPNT');
writeln ('Actual value of integral is',(fint(b)-fint(a)):7:4);
writeln ('n':6,'Approx. integral':29);
FOR i := 1 to nmax DO BEGIN
midpnt(a,b,s,i);
writeln (i:6,s:24:6)
END
END.